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Spin fluctuation-induced superconductivity in two-dimensional organic compounds such as k- 
(ET)2X is investigated by using a simplified dimer Hubbard model with right-angled isosceles 
triangular lattice (transfer matrices — t, — r'). The dynamical susceptiblity and the self-energy 
are calculated self-consistently within the fluctuation exchange approximation and the value for 
Tc as obtained by solving the linearized Eliashberg-type equations is in good agreement with 
experiment. The pairing symmetry is of A^2_y2 type. The calculated ((7/r)-dependence of Tc 
compares qualitatively well with the observed pressure dependence of Tc. Varying the value for 
r'/r from to 1 we interpolate between the square lattice and the regular triangular lattice and 
find firstly that values of Tc for k-(ET)2X and cuprates scale well and secondly that Tc tends 
to decrease with increasing r'/r and no superconductivity is found for r'/r = 1, the regular 
triangular lattice. 
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Major topics in condensed matter physics in recent 
years may certainly include the unconventional super- 
conductivity observed in high-Tc cuprates, heavy elec- 
tron systems and organic compounds. Among various 
possible mechanisms proposed so far the spin fluctua- 
tion mechanism seems to have possibilities of explaining 
all of these three. El The superconductivity in heavy elec- 
tron systems is believed to be almost certainly due to 
the spin fluctuation mechanism. Although the mecha- 
nism for cuprates is still controversial, the spin fluctua- 
tion mechanism seems to be rather unique in its ability of 
explaining not only the values of Tc and the pairing sym- 
metry but also various anomalous physical properties in 
the normal state as well as in the superconducting state 
in a systematic fashion at least in the optimal and over- 
doped concentration regimes;© the underdoped regime 
still remains to be clarified. 

As for the two-dimensional (2D) organic superconduc- 
tors the spin fluctuation mechanism seems to be the only 
available mechanism provided the superconducting or- 
der parameter is anisotropic, ^^iftf d-wave, as was indi- 
cated by recent investigationsooEP Major differences of 
this problem from that of high-Tc cuprates are that the 
superconductivity occurs without doping and in many 
cases in the metallic side of a metal-insulator Mott tran- 
sition and thus the system should be in the intermediate 
correlation regime. For example, the t-J model should 
safely be excluded. 

We wish to discuss the spin fluctuation mechanism of 
superconductivity in 2D organic compounds, keeping n- 
(ET)2X [ET = BEDT-TTF,X = Cu{N(CN)2}X',X' = 
CI, Br] in mindfl-' According to experiment k-(ET)2A" 
[X' ~ CI] is an antiferromagnetic insulator and under 
increasing pressure iJ. undergoes a insulator to super- 
conductor transitionE' Each layer of molecules in these 



compounds may be regarded as consisting of dimers each 
of which has one hole in the antibonding dimer orbital 
of the highest occupied molecular orbitals (HOMO).Bla' 
There are transfer matrices — r (r > 0) between dimers 
and the additional consideration of intra-dimer electron 
interaction U naturally leads to the Hubbard model. For 
U ^ T we have antiferromagnetic ground state and with 
decreasing U/t we encounter an insulator to metal Mott 
transition at a certain value of U/t. 

We now wish to study possible superconductivity on 
the metallic side of this transition. For this purpose 
we make use of a half-filled single band Hubbard model 
consisting of antibonding dimer orbitals with the inter- 
dimer transfers and the intra-dimer electron interaction. 
The spin fiuctuations are treated within the renormal- 
ized random phase approximation (RRRA) or the fiuctu- 
ation exchange (FLEX) approximation.EiJ Although the 
formalism of the self-consistent renormalization ISCR) 
theory takes into account the vertex correctionsJiJ^ as- 
sociated numerical jobs seem to be too heavy to carry 
them through. The results of FLEX for cuprates so far 
seem to be fairly successful and there are a few argu- 
ments Qn_ii.ertex CDrieclions in support of the FLEX ap- 
proachiBBBBBe 

Thus we see that the problem itself is reduced to the 
one quite similar to the Hubbard model description of 
cuprates. Only difference is the lattice structure and the 
transfer matrices. As a matter of fact we may simplify, 
to a good approximation, the model for k-{EiT)2X to a 
square lattice with the nearest neighbor transfers r and 
one of the cross diagonal second neighbor transfers r' 
(say up right corner to down left), thus making a right- 
angled isosceles triangular lattice with transfer matrices 
— r for the two sides and — r' for the base. Fig. l(a).ta) 
The value for r'/r as estimated from the presently ac- 
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Fig. 1. (a) The model unit cell and the transfer integrals, (b) 
Unperturbed Fermi surface for t'/t = 0.8. Dashed lines show 
the antifcrromagnetic zone boundary. 



cepted values for the. transfer integrals is about 0.8 (close 
to 1) for K-(ET)2XBEI' This reminds us of the magnetic 
frustration problem in an equilateral triangular lattice. 
So we study not only the case of k-{ET)2X but also the 
problem for varying values of t'/t, including the regular 
triangle: t'/t = 1. 

The model Hamiltonian is given by 



(q,if2,„) = [/ + [/2 



3 1 

-^Xs (q,ii^m) - ^Xc (q,i^2m) 



'2 {^s(9:i^?m) - Xc 



(9) 



and 



H 



Xs (g,i^m) 
Xc (g,ii^m) 

T 



Xs (g. i^r, 



Xc (g,i^m) 



(10) 



Xs (g, i^?m) = - — [G (fc + q, icc)„ + ii7„) G (fc, itj„) 

iVp — ^ 



(1) 

(2) 



c k 



with 



-2r (cos fcajO + cos kyo) 

—It' cos (fc^a + hyo) 



(3) 



+i^(fc + q,iw„ + iJ7„0i^(fc,it^n)] , (11) 
Xc (g, i^^m) = -TT [G (fc + q, ici^n + ii7m) G (fc, itJ„) 

-F {k + q, iujn + ifim) F {k, itj„)] , (12) 

where w„ = (2n + 1) ttT and Slrn = 2rmrT are the Fermi 
and the Bose Matsubara frequencies, respectively and iVo 
is the number of dimers in the crystal. We study have 



where (ij) and (ij) indicate the nearest neighbor pahs the possibility of singlet pahing. 



and the second neighbor side diagonal pairs, respectively, 
as shown in Fig. 1(a) and a is the lattice constant. 

The equations for the normal and anomalous Green's 
functions are given by 

G(fc,iw„) = G(") (fc,iw„) 

+G(°) (fe, iujn) [s^^^ (fc, iujn) G (fc, iLUn) 

~S^^'> {k,icun)F* (fc,iw„)] , 
F* (fc,ia;„) = G(")(-fe,-icc>„) 

X (-fc,-iw„)F*(fc,iLj„) 

+S'^^A~k,-iij^)G{k,iu;n)], (5) 

where G^°-' (fc,ia;„) is the Green's function for a non- 
interacting system and the self-energies due to the spin 
and charge fluctuations are given as follows: 

^ .^1^(1) (g,ir2„) 



Confining ourselves here to the transition temperature 
Tc and the normal state properties we may linearize the 
equations (4, 5) and (11, 12) with respect to the anoma- 
lous self-energy S^^^ (A;,iaj„) or f (fc,i[j„). Tc can be 
calculated as the highest temperature where the follow- 
ing equation for the normalized anomalous self-energy 
/ (fc, iuJn) = ^^^^ {k, iujn) |G (fc, icJn)] has a non-trivial so- 
lution: 



(4) /(fc,icc>„) 



= -^EiG(fc,-n)i 



1 



X |G(p,icJm)| /(p,iw„) . 



(13) 



Z-fl) {k,iL 



Nn 



xG(fc - g,iw„ - ii7m) , (6) 



The symmetry of the order parameter is given by that of 
/ (fe, iujn)- Since the kernel of eq. (13) has a full symme- 
try of the model Hamiltonian we may classify the order 
parameter according to the irreducible representations 
of the symmetry group. There are four one-dimensional 
irreducible representations Ai, A2, Bi and B2. We take 
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Fig. 2. Density of states for various values of U/t. 



Fig. 5. Tc vs. U/t. 
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Fig. 3. Quasi-particle dispersion. 
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Fig. 4. Dynamical susceptibility. 



advantasa,of using this symmetry property in the cal- 
culation.ll2P For numerical calculations we take 64 x 64 
fc-points and 512 and 2048 Matsubara frequencies for 
susceptibility and self-energy, respectively. 

We first show the results of calculation for k-(ET)2X 
{t' /t = 0.8). The unperturbed Fermi surface is shown 
in Fig. 1(b). The one particle density of states and 
the quasi-particle dispersions are shown in Figs. 2 and 



3, respectively, together with the corresponding values 
for U = 0. The dynamical susceptibilities are shown in 
Fig. 4 for various frequencies. We see that incommensu- 
rate peaks around (tt, tt) and an asymmetric peak around 
(tt, — tt) are strongly enhanced. 

The superconducting transition temperature Tc is 
given as Tc — O.OlOSr for [/ = 9t and the order pa- 
rameter has the A2 symmetry {x^ — y^) . Estimating 
r ~ 0.07eV, we get Tc ~ 9K in fair agreement with 
the observed value ~ lOK. Fig. 5 shows the plot of 
Tc vs. U/t. Tc shows a weak maximum aX U /t % and 
then decreases with decreasing U/t. Since U/t should 
decrease with increasing pressure, this tendency is con- 
sistent with the observed pressure dependence of Tc if we 
assume a proper critical value {U /t)^, sauaxound 10, for 
the insulator-superconductor transition.LrEJ According 
to the mean field calculation, which normally underesti- 
mate the critical point for transition, the first order tran- 
sition between the ip^^ulator and metallic phases takes 
place at U/t = A.2n It still remains to calculate the 
critical value (U/t)^ by using the same FLEX approxi- 
mation for both antiferromagnetic and superconducting 
states at T = 0. 

We also note that the calculated values of Tc for k- 
(ET)2X and cuprates properly scale. A crude model 
for some of the cuprates, say LSCO, may correspond 
to the case of r'/r m with less than half-filled car- 
riers. According to the previous calculations we have 
Tc/r = 0.0273 (0.021) for t'/t = (0.15 for cross di- 
agonal traBsfers), U/t = 6(4) and the carrier number 
n = 0.875.li3'li3) Since the band width of cuprates is con- 
sidered to be roughly 2 ~ 4 times larger than that in 
K-(ET)2X, the calculated relative values of Tc compare 
rather well with experiment. 

Next we discuss the (r'/T)-dependence of Tc for the 
half-filled case. It is interesting to find that eq. (13) 
has no solution for t'/t — 1 (the regular triangular lat- 
tice) in any reasonable range of the value for U/t, or 
U/t < 16, the highest value studied. For U/t » 1 and 
t'/t = 1 we have an antiferromagnetic Heisenberg model 
with a regular triangular lattice, a famous frustrated sys- 
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Fig. 6. Wave vector dependence of the susceptibility for t'/t 
1.0, the regular triangular lattice, for various values of U/t. 



teni.0 Around the Mott transition under pressure the 
local moments disappear and the metallic phase seems 
to be characterized by a wave vector dependent mag- 
netic susceptibility with three broad peaks as is shown 
in Fig. 6. This situation does not seem favorable for 
the spin fluctuation-induced superconductivity as may 
be seen from a weak coupling argument Ej^ 

In the case of t'/t = 0.4, the antiferromagnetic peak 
of the dynamical susceptibility is more significant than 
for t'/t ~ 0.8 and the superconductivity appears for 
smaller U and the value of is higher as may natu- 
rally be expected. For example, we have T^/t — 0.014 
and 0.016 for U/t = 3.3 and 3.7, respectively. For a 
half-filled band we may conjecture as follows: Consider- 
ing a phase diagram in U/t vs. t'/t plane, the critical 
boundary for antiferromagnetism (f//T)^p is lower than 
that for superconductivity {U/t)^^ for small values of 
t'/t where the nesting condition is well satisfied. With 
increasing t'/t the frustration increases and its destruc- 
tive influence is stronger for antiferromagnetism than for 
superconductivity and thus {U/t)^^-, becomes lower than 
(J7/t)^p beyond a certain value of t'/t. For a definite 
conclusion it is necessary to compare the free energies of 
both of these states. 

As for the normal state properties the uniform sus- 
ceptibility is almost constant and has weak tendency of 
decreasing with decreasing temperature. The calculated 
values for J2q Inix (q, ^) /^, the quantity proportional to 
l/TiT, Ti being the nuclear spin-lattice relaxation rate, 
has a Curie- Weiss like temperature dependence down to 
Tc- On the other hand the experimental results on k- 
{ET)2X [X' = Br] show a substantial decrease of the 
susceptibility as Tc is appmajcied and the peak structure 
of 1/TiT much above Tc.BoEI' These behaviors remain 
to be explained. 

In summary we have studied the spin fluctuation 
mechanism of superconductivity in k-(ET)2X by using 
a FLEX approximation for a dimer Hubbard model with 
the effective transfer integrals — t, — t' deduced from 
the presently accepted values for the transfer integrals 
and varying values of U. We found superconductivity 
of dj.2_j,2 types. The value for Tc and its U/t (pres- 



sure) dependence compare well with experiments. We 
point out that if we vary the value of t'/t in the present 
model for k-(ET)2X between and 1 the square lat- 
tice and the regular triangular lattice Hubbard models 
are interpolated. We flnd that the values of Tc for k- 
{ET)2X and for cuprates scale fairly well and there is 
no superconductivity of singlet pairing for the triangu- 
lar lattice within a reasonable range of the value for U. 
However, the anomalous normal state properties in the 
uniform susceptibility and the nuclear spin-lattice relax- 
ation rate, bearing resemblance to the pseudo-spingap 
phenomena in cuprates, remain to be explained. Also, 
the transition between the antiferromagnetic insulator 
and superconducting phases is still to be investigated. 

We would like to thank Prof. K. Kanoda, Dr. S. 
Nakamura and in particular Dr. T. Takimoto for helpful 
discussions. 
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